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Abstract. We report a theoretical study of the quasiparticle lifetime and the quasiparticle mean free path 
caused by inelastic electron-electron scattering in ferromagnetic compounds of the Co-Fe and Ni-Fe systems. 
The study is based on spin-polarized calculations, which are performed within the G'^W° approximation 
for equiatomic and Co(Ni)-rich compounds, as well as for their constituents. We mainly focus on the spin 
asymmetry of the quasiparticle properties, which leads to the spin-filtering effect experimentally observed 
in spin-dependent transport of hot electrons and holes in the systems under study. By comparing with 
available experimental data on the attenuation length, we estimate the contribution of the inelastic mean 
free path to this length. 
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1 Introduction 

Spin-dependent transport of hot electrons and holes in 
ferromagnetic materials is one of the basic phenomena 
exploited in constructing spintronic devices. As an ex- 
ample, one can mention spin-valve and magnetic tunnel 
transistors, which are also used to experimentally study 
the spin- dependent transport in different many-electron 
systems pQ[51[31|3] . Such a study aims to solve the problem 
related with controllable search for desirable parameters of 
constructed devices. From theoretical side, it is important 
to determine a role played by inelastic electron-electron 
(e — e) scattering in formation of operating characteristics 
of the devices. To this end, the finite lifetime caused by 
the e — e scattering and the corresponding mean free path 
of electrons and holes should be analyzed as a function of 
exciting energy. 

For pure ferromagnetic metals, a similar analysis of the 
mentioned quasiparticle properties was made within both 
semiempirical approaches and ab initio calculations. Some 
of the calculations were performed within the G^W^ ap- 
proximation (for a recent review, see, e.g., Ref. [S]). This 
approximation is a non-self-consistent variant of the GW 
approach based on Kohn-Sham states evaluated within the 
local spin density approximation (LSDA) with or with- 
out a Hubbard U correction as a starting point 0. Also, 



the calculations were done within the LSDA+DMFT ap- 
proach [511^110] ■ which combines the LSDA and the dy- 
namical mean- field theory (DMFT). This approach con- 
tains two parameters (the averaged screened Coulomb in- 
teraction U and the exchange interaction J) , which are de- 
fined outside the scope of the approach and can be taken 
from experiment or calculated, e.g., within the constrained 
local density approximation or the constrained random 
phase approximation (see, e.g., [TT1[T^[TB] V In that sense, 
the Giy-|-DMFT approach [T3] can be considered as a pa- 
rameter free one. 

However, in the most cases ferromagnetic alloys Coa;Fei_ 
and NijjFei-a; with large x are used in spintronic devices. 
To all appearance, the first attempt to theoretically study 
the lifetime and the inelastic mean free path in ferromag- 
netic alloys and compounds has been made by the authors 
in [15[ . In the work cited, the properties of quasipartil- 
ces (electrons and holes) in the ordered (B2) and disor- 
dered (body-centered cubic) CoFe have been calculated 
from first principles within the mentioned G^W'^ approxi- 
mation. Also, the contribution of the localized d states to 
the above properties has been analyzed by considering the 
Coulomb interaction screened by only d electrons instead 
of the fully screened Coulomb interaction. 

In the present paper, we proceed with the theoretical 
study of the properties of quasiparticles in ferromagnetic 
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compounds, as well as in their constituents. In order to 
be closer to the Co- and Ni-rich ferromagnetic alloys used 
in practice, among the compounds we consider CosFe and 
NiaFe with the DO3 and LI2 structure, respectively. We 
perform ab initio calculations within the G^W^ approxi- 
mation with the fully screened Coulomb interaction found 
in the random phase approximation (RPA). To demon- 
strate possible effects that structure and stoichiomctry 
changes may have on the quasiparticle properties under 
study, we additionally consider the ferromagnetic com- 
pound NiPe with the LIq structure. For completeness of 
the picture of such an effect, we also analyze the results 
obtained in [15] for the B2 CoFe. 

The paper is organized as follows. In Sec. [H we briefly 
describe the scheme that is used to evaluate the spin- 
resolved quasiparticle lifetime and inelastic mean free path 
within the G'^W^ approximation. Also, we put here some 
calculation details. In Sec. [31 we present our main results of 
extensive calculations carried out for quasiparticle proper- 
ties in the mentioned compounds and pure metals. On the 
base of these results, we analyze how changes in structure 
and stoichiometry modify quasiparticle properties evalu- 
ated from the G'^W'^ calculations. Finally, the conclusions 
are given in Sec. |4| 



2 Approximations and calculation details 

In this section, to make the paper self-sustained we briefly 
describe the G°W° approximation [T^ to the quasiparti- 
cle self-energy, which is used in the paper. Unless stated 
otherwise, atomic units are used throughout, i.e., e^ = 
h ~ m ^ 1. 

The G'^W'^ self-energy of a quasiparticle with the spin 
a is deflned as 

i:,(ri,r2;c^) ^ ^ J du'e"^^' G'i{ri,r2;uj') 

X M^0(ri,r2;w-c^'), (1) 

where the convergence factor exp(ir/aj') indicates that the 
integration contour is closed in the upper half-plane of a;'. 
The Green function entering Eq. ([!]) is deflned as [T7] 



G0(ri,r2;c.)-^ 



V'k«a(ri)7/'^„^(r2) 
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and corresponds to the Kohn-Sham equation with the 
exchange-correlation potential V^'~^ obtained within the 
LSDA. In Eq. ([5]), the positive inflnitesimal S characterizes 
the way of going around poles in the course of integra- 
tion. The states {i^knaj^Wna} are calculated by the self- 
consistent tight-binding linear muffln-tin orbital method 
[T5] with the use of the crystal potential constructed within 
the atomic-sphere approximation. 

The screened Coulomb interaction W'^ participating 
in the definition of the self-energy ([!]) is found within the 



RPA, where the irreducible polarizability P'^ is expressed 
as 

P0(ri,r2;c.) = ^ ^^ y'dc.'GO (n, r2;c.') 

xG^(ri,r2;w'-fc.). (3) 

Owing to the translation symmetry, the irreducible polar- 
izability can be expanded into a series [TTlfT^ 

P0(ri,r2;c^) =^Bk.,(ri)/^'5-(k,a;)i?£^.(r2), (4) 

where {i?ki} is a set of basis functions, which satisfy the 
Bloch theorem and are normalized to unity in the unit 
cell volume f2. These basis functions arc constructed us- 
ing pair products of linear muffin-tin orbitals localized at 
the same lattice site [1D|. The irreducible polarizability 
Pjj(k, a;) is derived with the use of the spectral function 
representation. The corresponding spectral function 

^iji^j'^) = /_^ 2^ (/k+qn'o- — fkncr) (BqjV'kno- 1 V'k+qn'cr) 
kcr nn' 

where /kno- is the Fermi factor, is obtained by substituting 
the Green function of Eq. ^ into the expression ((3|) for 
P" and using the expansion ^ along with the relation 
ImP°(k, w) = — 7rS'° (k, w)sgn(w). In calculations, the S- 
function is replaced by the Gaussian exp{—uj'^/j^)/{j^/n) 
with 7 = 0.136 eV [21]. 

The inverse quasiparticle lifetime or decay rate is de- 
fined as T^^l^ = 2Zk„o-|Inirk„<T(ekr 
malization factor Zkno- is given by 



where the renor- 



■^knfT — 



1- 
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With the use of the expressions presented above, the imag- 
inary part of the self-energy matrix elements Skna-i^^) = 
('0kncr|^(T(w)|'0kncr) entering the definition of the quasi- 
particle decay rate takes the form 

ImSl^nai^) = Ty^y^(V'k»<T^q-k»'<T|-Bqi:) 
qn' ij 

X ImM^° (q, ±eq_kn'<T T w) (6) 

X (-Bqj|V'q-knvV'kncr)G'(±eq-k«'(T T w), 

where the upper (lower) sign corresponds to the excit- 
ing energy w < Ep (cj > Ep) and the sum over occu- 
pied (unoccupied) states. Here, Ep is the Fermi energy. 
In Eq. ([5]), 0{x) is the step function and Wij are matrix 
elements of the screened interaction in the {Pki} basis, 
which are defined by Pj^„(k, w) and Coulomb-interaction 
matrix elements [17] . Having obtained the imaginary part 
of Z'kncr('i^); the real part of the latter is found from the 
Hilbert transform. 

Both in Eq. ^ and in Eq. ([5]), the s, p, and d bands are 
involved in the sums over the band indices. As to the sums 
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over momenta k, a set of points 20 x 20 x 20 is used for the 
face- and body-centered cubic (fee and bee) structures and 
the B2, LI2, and DO3 structures. For the LIq structure 
with the primitive vectors (-1/2, 1/2, 0), (1/2, 1/2, 0), and 
(0, 0, 1), we use a set of points 20 x 20 x 12. The quasipar- 
ticle hfetimes are calculated at all k-points of these sets, 
which belong to the irreducible Brillouin zone, what en- 
sures the momentum averaging of the lifetimes to be well 
converged. The number of optimal product basis functions 
per atom is 40. 

As regards the well-known problem with the slow con- 
vergence of ReZ'kno- with respect to the number of unoc- 
cupied states and basis functions, it is worth noting that 
in order to consider G'^W'^ corrections to the LSDA band 
structure, a minimal requirement in our case is to involve / 
band and to increase the number of optimal product basis 
functions per atom to 115 with the inclusion of core states 
(see, e.g., [TTll^ ). Since we are interested in the quasipar- 
ticle lifetime determined by the self-energy, where W'^ is 
computed within the RPA instead of, e.g., the plasmon- 
pole model, the mentioned modifications of the calculation 
parameters make a study of the lifetime in the ferromag- 
netic compounds hardly realized. At that, as an analysis 
has shown by the example of Ni, ImZ'tncr and Z^na con- 
verges noticeably faster than ReZ'kna- At the calculation 
parameters used in the paper, the imaginary part of the 
self-energy and the renormalization factor do not differ 
substantially (within 4 — 7% on average in the considered 
energy range) from these obtained within the requirement 
just mentioned. At least, the found differences do not af- 
fect the issues, which we discuss below. 

In order to estimate the inelastic mean free path (IMFP) 
of quasiparticles as A^~^ = Ivkno-lTkna, the quasiparticle 
velocity vit^na can be evaluated from the expression (see, 
e.g., Ref. [IS]) 



Vkr: 



■VkReziZ'k„<x(a;)|. 



(7) 



where v^^^ = Vkfikncr and the difference Z\Z'kncr(w) = 
Sknai^^) - {ipkncr\Vj^'-^\iJkna), which is treated as an ex- 
pansion parameter. In the present calculations, for sim- 
plicity we neglect the second term in the right hand side 
of Eq. (O. This leads to the IMFP given by the frequently 
used formula 



a: 



kncT 



2ImJ7kno-(ekncr) 



(8) 



It is worth noting that this formula, which, e.g., in |23j 
has allowed to obtain theoretical results in good agree- 
ment with experimental data in the case of beryllium, 
does not contain the real part of the self-energy matrix 
elements. On the one hand, this means that many-body 
effects cannot be taken into account in full measure. But 
on the other hand, such an IMFP is free of the conver- 
gence problem mentioned above and is completely defined 
by the underlying (LSDA, in our case) band structure. For 
a pragmatic point of view, the used approximation seems 
to be a quite accurate and feasible method for predicting 
the quasiparticle lifetime and the IMFP in ferromagnetic 
metals and compounds on the same footings. 
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Fig. 1. Density of states in the ferromagnetic compounds CoFe 
and CoaFe with the B2 and DO3 structure, respectively, as well 
as in the bcc Fe and the bcc Co with the lattice parameter of 
CoFe. 



3 Results and discussion 

In the present calculations, the cubic lattice parameter a 
for CoFe (B2) and CoaFe (DO3) is chosen to be equal to 
2.856 Agll and 5.657 A[15], respectively. The same lattice 
parameter as for CoFe is used for the bcc Fe and the bcc 
Co. This parameter is very close both to the equilibrium 
parameter of pure iron and to the parameter of bcc cobalt 
films (2.82 ± 0.01 Agn]) grown on GaAs. Moreover, the 
parameter a of CosFe with the bee-based DO3 structure 
is about doubled a of the mentioned bcc Co films. In NiFe 
(Llo), the parameter a has the value of 3.585 A, which 
was estimated by linear interpolation between 57.0 and 
44.3 at. % Ni in Nii_^Fe^ fee alloys [24]. For NigFe with 
the LI2 structure and for the fee Ni, the lattice parameter 
is taken to be equal to 3.552 A|24j. 



3.1 Co-Fe system 

Fig-fflshows the calculated density of states (DOS) both in 
the considered ferromagnetic compounds of the Co-Fe sys- 
tem and in the pure metals forming these compounds. The 
presented curves demonstrate the behavior of the DOS 
as a function of energy, which corresponds to the LSDA 
band structure. We plot these curves with the purpose to 
explore specific features, which already at the stage of an- 
alyzing the DOS allow one to qualitatively describe the 
behavior of the ratio t-^/t^ between the lifetimes of quasi- 
particles with different spin orientations and how this ra- 
tio changes as we move from Fe, to CoFe and CoaFe, and 
then to Co. Such a description is based on the fact that 
the DOS directly depends on dispersion of ekn that along 
with ip'kn determines quasiparticle lifetimes. 

Due to the bcc-based structure of the considered ma- 
terials, the DOS has a form that is typical for the bcc 
transition metals: a profound minimum halves the d-band 
[27] . Moreover, the DOS in the bcc Fe is similar to that 
in the bcc Co, where the main effect caused by moving 
from iron to cobalt is a change in Ep due to an additional 
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Table 1. Calculated local magnetic moments (in hb) 



Material 


/iPc 


MCo 




Co (bcc) 




1.72 


This work 






1.3-1.7 


Exp. mmm 


CoaFe (DO3) 


2.62 


1.72 / 1.71 


This work 




2.61 


1.78 / 1.78 


Theory |34| 


CoFe (B2) 


2.73 


1.74 


This work 




2.82 


1.74 


Theory [35] 




2.92±0.02 


1.62±0.02 


Exp. [36] 


Fe (bcc) 


2.21 




This work 




2.22 




Exp. [37] 



electron. As a consequence, the densities of states in Co, 
CoFe, and CosFe are very close to each other, especially in 
the spin-up subsystem. As to the spin-down subsystem, a 
distinguishing feature of these cobalt-containing materials 
is the DOS shape in the vicinity (approximately ±1 eV) 
of the Fermi level. In this sense, an important difference 
between, e.g., CosFe and CoFe is that the former has the 
DOS (predominantly formed by d states) that is higher 
than that in the latter. This means that in the vicinity 
of Ep the DO3 CosFe is more spin polarized than the B2 
CoFe. 

The obtained LSDA band structure is characterized 
by the local magnetic moments, which are listed for com- 
parison with experimental and other theoretical data in 
Table [T] Note that in the case of the DO3 CosFe two 
values of /^co are presented. It caused by the fact that 
there are two inequivalent sites for Co atoms in the DO3 
structure. As is seen from the table, in the considered 
compounds the magnetic moment for Co remains prac- 
tically unchanged as compared with that in the pure bcc 
cobalt. At the same time, the magnetic moment for Fe be- 
comes noticeably larger upon moving from the pure bcc 
iron to the compounds. Such a behavior of fico and ^fc is 
in agreement with available experimental data. 

As noted above, the mentioned changes in filling, width, 
and shape of d-band should affect the lifetime as a func- 
tion of exciting energy. In the case of paramagnetic simple, 
noble, and transition metals, such an influence has been 
considered in [2811291(30] in details. In the present work, 
we can additionally analyze the ratio t-^/t^. Actually, as 
is seen from Fig. [TJ in the vicinity of Ep (especially for 
positive exciting energies) the density of spin-up states is 
considerably smaller than the density of spin-down states. 
As a result, we can expect that, at least in this energy 
rage, the ratio t-\-/ti will be substantially greater than 
unity. An exception is constituted by the bcc iron that 
should be characterized by the above ratio being smaller 
than unity. 

The results of our calculations of the quasiparticle prop- 
erties in the B2 CoFe are shown in Fig. [5] In this figure, we 
demonstrate both momentum-resolved and momentum- 
averaged lifetime and renormalization factor for spin-up 
and spin-down quasiparticles. In order to clearly repre- 
sent the exciting-energy dependence of the quasiparticle 
lifetime, the so-called scaled lifetime (that is the lifetime 




Fig. 2. Scaled lifetime rk„o- x (ek^o- — Ep)^ (upper panel) and 
renormalization factor Zj^n^ (lower panel) as functions of the 
exciting energy ekntr — Ep for spin-up (light gray points) and 
spin-down (dark gray points) quasiparticles in the B2 CoFe 
ferromagnetic compound. Solid and dashed lines show the cor- 
responding scaled lifetimes averaged over momentum k for a 
given exciting energy. 




Fig. 3. Same as in Fig. [2] but in the DO3 CoaFe ferromagnetic 
compound. 



multiplied by squared exciting energy) is shown. Analyz- 
ing the data presented in Fig. [3] we would like to note 
the significant difference between the lifetime of electrons 
with spin up and spin down. Such a difference is one of 
the main reason of appearing the so-called spin-filtering 
effect that is observed in spin-dependent transport of elec- 
trons in ferromagnetic materials |31[S] . The lifetime of holes 
(especially with spin down) depends strongly on the mo- 
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mcntum k at a given exciting energy. It follows from the 
observed spread in values of Ti^na at euna < Ep . 

As regards the renormalization factor characterizing 
the spectral weight of quasiparticles, like the densities 
of states for corresponding spin subsystems the factors 
Zknt and Zknj. as functions of exciting energy are quite 
similar, in many respects, but with some shift on energy 
scale, which is determined by the exchange splitting of the 
bands. For both spins, there is a "depression" in the en- 
ergy range that corresponds to the dominant contribution 
of the localized d states to the total DOS. This means 
that many-body effects have the most impact on these 
states (see, e.g., Ref. |3H]). In the energy range from ^0 to 
~3 eV, where a comparatively high density of spin-down 
states corresponds to a low density of spin-up states, one 
observes the largest difference between Zi^„-^ and Zi^ni- 
At higher energies, the renormalization factor of electrons 
depends weakly on spin. 

Fig. [3] shows our results on the quasiparticle lifetime 
and renormalization factor in the case of the ferromagnetic 
compound CoaFe with the DO3 structure. As a whole, 
for electrons the scaled lifetime (momentum-averaged and 
momentum-resolved) differs slightly from that in the B2 
CoFe, whereas for holes in the exciting-energy range from 
to -3 eV the quantity under consideration demonstrates 
noticeably smaller values and substantially weaker depen- 
dence on the momentum k at a given exciting energy. 
At that, in the vicinity of the Fermi energy, the scaled 
lifetime of spin-up holes is longer than that of spin-down 
holes practically to the same extent as the scaled lifetime 
of spin-up and spin-down electrons. 

Upon moving from the B2 CoFe to the DO3 CoaFe, the 
renormalization factor undergoes a visible change in mo- 
mentum dependence at a given energy (especially at Ep) 
and, being momentum-averaged, in dependence on excit- 
ing energy. At that, the mentioned correlation between 
values of the renormalization factor and the d-states con- 
tribution to the DOS remains evident as before. Note that 
in this case Zk„i has a profound minimum at the Fermi 
level. This minimum is caused by a quite narrow and high 
peak located practically at Ep in the density of d states. 
The factor Zknt demonstrates a similar minimum at about 
— 1.5 eV, which is also correlated with a peak in the den- 
sity of d states (see Fig. [1]). 

In order to demonstrate what effect the changes in 
structure and stoichiomctry have on the lifetime of quasi- 
particle in the considered ferromagnetic materials, we show 
all the obtained results on the Co-Fe system in Fig. 2] As 
is evident from the figure (see the upper panel), the life- 
time of spin-up electrons in the compounds is longer than 
that in their constituents. Spin-down electrons have the 
longest lifetime in the bcc iron practically in the whole 
considered energy range. For holes, the longest lifetimes 
are observed in the B2 CoFe in both spin subsystems. 
Note the closeness of the hole lifetimes in the DO3 CoaFe 
and in the bcc cobalt. In fact, in all the cases, we deal 
with a spin asymmetry of the quasiparticle properties, 
which can be represented by the already mentioned ra- 
tio of the lifetimes of spin-up and spin-down quasiparti- 
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Fig. 4. Scaled lifetime averaged over momentum k for a given 
exciting energy (upper panel) and ratio r-f/ri between the aver- 
aged lifetimes of quasiparticles with different spin orientations 
(lower panel) as functions of exciting energy for the ferromag- 
netic compounds CoFe and CoaFe with the B2 and DO3 struc- 
tures, respectively, as well as for the bcc Fe and the bcc Co 
with the lattice parameter of CoFe. Inset in the lower panel 
shows the inverse ratio r^/r-|- for holes. 
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Fig. 5. Dependence of the momentum-averaged IMFP A|^~^ 
on exciting energy. The experimental data on the attenuation 
length for spin-up and spin-down electrons in Co84Fei6 are 
taken from [3] 



cles. A deviation of the ratio from unity reveals the spin 
asymmetry of the considered quantity. In Fig. SI we also 
plot the ratio t-^/ti between the spin- up and spin-down 
momentum- averaged lifetimes (see the lower panel). One 
can see that the expected values qualitatively estimated 
from the shape of the DOS in the vicinity of the Fermi 
level are confirmed by the results presented in the figure. 
Actually, in the exciting-energy range from to ~0.7 eV 
the ratio t^/t^ amounts to 8 for the B2 CoFe and ex- 
ceeds 12 for the DO3 CosFe and the bcc cobalt. For the 
bcc iron, the ratio varies from ^--^0.2 to ^0.8 in the above 
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energy range and only starting from ~^1.2 cV becomes 
greater than unity. However, for holes in Fe the inverse 
ratio Ti/T-|- depends weakly on energy and comes to '^2 
(see the inset in the lower panel of Fig. 2]). It is worth not- 
ing that in contrast to the bcc iron and the bcc cobalt the 
considered ferromagnetic compounds are characterized by 
the ratio t^/t^, which over the wide energy range (from 
~ 1 to 5 eV) is large (about six) and varies slightly with 
exciting energy. For electrons with energies less than ~ 1 
eV, the DO3 CosFe provides the largest t^/t]^. 

As regards the IMFP expressed as the product of the 
lifetime and the quasiparticle velocity, an analysis has 
shown (see, e.g., Rcf. [15]) that a noticeable difference 
between spin-up-electron and spin-down-electron veloci- 
ties increases the revealed significant difference between 
the lifetime of electrons with spin up and spin down. This 
means that for the ferromagnetic materials under study 
the calculations within the G'^W^ approximation predict 
a strong spin- filtering effect, which manifests itself in a 
giant spin asymmetry of the IMFP as is clearly seen from 
Fig. [S] However, the experimental data (also presented in 
Fig. [S]) on the spin-dependent transport of electrons in 
Co84Fei6 films on the GaAs(OOl) surface |3] indicate that 
the mentioned effect is not so strong. At that, the calcu- 
lated momentum-averaged IMFP for spin-down electrons 
in the compounds and in the bcc cobalt is quite close 
to the corresponding attenuation length and nicely repro- 
duces the weak energy-dependence of the latter, while the 
difference between the experimentally observed attenua- 
tion length and the calculated IMFP for spin-up electrons 
is substantial. 

To analyze the discrepancy between the theoretical 
and experimental data, first, we should note that the IMFP 
calculated within the G'^W^ approximation does not in- 
clude the contribution of the decay channels caused by 
spin fluctuations. As was shown in J39[|40[HT] for param- 
agnetic and ferromagnetic materials, taking these fluctu- 
ations into consideration leads to a decrease in the quasi- 
particle lifetime and, as a consequence, in the IMFP \%~'^ . 
At that, AS"*^ appears to be most sensitive (see, e.g., Ref. [5^ 
Next, following [5] we represent the measured attenua- 
tion length as a result of the sum of two terms: 1/A^^'' = 
1/A^~*^-|-1/A|^*'"°; where, in addition to the inelastic electron- 
electron scattering contribution (the calculated A^"*^), there 
is a contribution that includes the terms coming from 
quasiclastic scattering by phonons and spin waves and 
from elastic electron scattering by defects and impurities. 
The latter depends on type of defects and their number 
and represents the term, of which contribution to the spin 
asymmetry of the attenuation length is expected to be 
minimal (see, e.g., Ref. [12])- The contribution of electron- 
magnon scattering is known to be significant for spin-down 
electrons (see, e.g., [13]). The term caused by electron- 
phonon scattering, has an effect on both spin-down and 
spin-up electrons. The dependence of this term on excit- 
ing energy is almost completely determined by the elec- 
tron velocity. This means that the corresponding decay 
rate is nearly constant. To illustrate how the inclusion of 
a contribution caused by a decay process weakly depen- 
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Fig. 6. Density of states in the ferromagnetic compounds NiFe 
and NigFe with the Llo and LI2 structures, respectively, as well 
as in the bcc Fe and the fee Ni. 



dent on exciting energy can affect the attenuation length 
of spin-up electrons, we approximate 1/AS^*''° as F^/v^ 
with r^ = 50 meV. In the order of magnitude, this F^ 
is close, e.g., to the electron-phonon broadening observed 
for spectral lines in different metals (see, e.g., I44]). The 
obtained results are presented in Fig. [SJ In the case of 
the bcc cobalt, we have the least difference between the 
experimental and theoretical data. To this, the behavior 
of the attenuation length as a function of exciting energy 
is nicely reproduced. Thus, in order to reach a satisfac- 
tory agreement with the experimental data, at least the 
electron-phonon broadening should be included. 



3.2 Ni-Fe system 

In Fig. [HI we show the results of the DOS calculations for 
two compounds of the Ni-Fe system and the pure fee nickel 
within the LSDA. The fee nickel and the compounds have 
the density of spin-up states with the continuous d-band 
)that is peculiar for fcc-based structures. At that the Fermi 
level is situated quite close to the edge of the d-band. 
The density of spin-down states possesses the fcc-typc ri- 
band in the case of fee nickel only, where the Fermi level 
corresponding to a peak at the edge of the d-band makes 
this ferromagnetic metal to be highly spin polarized. The 
compounds NiFe and NisFe demonstrate the DOS with 
the minimum at Ep, which is rather close to the ordinary 
half- filled bcc d-band (especially in the LIq NiFe case). 
This means that at least in the vicinity of Ep the fee nickel 
should demonstrate a quite large ratio t^/t^ as compared 
with the compounds. 

The found local magnetic moments are presented in 
Table [D As is evident from the table, again, as in the 
case of the Co-Fe system, in compounds the magnetic mo- 
ment for Fe becomes larger than in pure iron. At that, the 
largest moment ^po is observed in the LI2 NiaFe. As to Ni 
atoms, with good agreement with experimental data the 
magnetic moment for Ni remains unchanged upon moving 
from the fee nickel to the compounds. 



I. A. Nechaev, E. V. Chulkov: Quasiparticle dynamics in ferromagnetic compounds of the Co-Fe and Ni-Fe systems 



Table 2. Calculated local magnetic moments (in ^b) 



Material 


/iPc 


MNi 




Ni (fee) 




0.61 


This work 






0.62 


Exp. [37] 


NigFe (LI2) 


2.85 


0.62 


This work 




2.87 


0.62 


Theory 1451 




2.97±0.15 


0.62±0.05 


Exp. m 


NiFe (Llo) 


2.59 


0.61 


This work 




2.6 


0.6 


Theory [47] 



Now, we turn to the analysis of the lifetime and the 
renormalization factor in the ferromagnetic compounds of 
the Ni-Fe system. In Fig. [7] we demonstrate both momentum- 
resolved and momentum-averaged lifetime (or, more pre- 
cisely, the scaled one) and the renormalization factor for 
spin-up and spin-down quasiparticles in the Llo NiFe. The 
main feature we would like to point out is a spread in val- 
ues of the scaled lifetime at a given exciting energy in the 
energy rage ±2 eV, what we did not observe in the Co-Fe 
system. Also, it is easily seen that the scaled lifetime of 
spin-down quasiparticles is quite symmetric with respect 
to zero exciting energy. In addition to the spread of values, 
the difference between the lifetimes of spin-up and spin- 
down electrons is not so large as, e.g., in the B2 CoFe. 
The same one can say about spin-up and spin-down holes. 
In the energy range from —1 to —3.5 eV, the lifetime of 
spin-up holes is very close to that of spin-down holes. 

As regards the renormalization factor, for spin-up quasi- 
particles this quantity averaged over k is closely approxi- 
mated to that in the B2 CoFe, apart, maybe, from more 
smooth behavior as a function of exciting energy due to 
the continuous d-band. For spin-down quasiparticles, the 
momentum-averaged Zknj, has two depressions (symmet- 
ric with respect to the Fermi level) owing to the profound 
minimum in the d-band at Ep, which separates predom- 
inantly iron states (above Ep) and predominantly nickel 
states (under Ep). This separation means that in the spin- 
down subsystem for electrons the probability to be in iron 
d states exceeds that to be in nickel d states, while for 
holes - vice verse. 

The obtained results on Tk„a and Zi^na in the LI2 
NisFe are shown in Fig. \7\ The presented data denote 
that as the nickel content increases the spread in values 
of Tkno- at a given exciting energy and different momenta 
tends to be narrowed. Around the Fermi level, on average 
the lifetime of spin-up and spin-down quasiparticles be- 
comes shorter. In addition to this, for holes in the energy 
range from —1.5 to —3.5 the momentum- averaged lifetime 
does not show any noticeable dependence on spin. For 
electrons, the difference between the momentum-averaged 
Tknt and Tknj, is rather close to that in the compounds of 
the Co-Fe system than to the difference in the Llg NiFe. 

As can be expected from the DOS shape, the renor- 
malization factor Zkn-|- does not undergo visible changes 
upon moving from the Llg NiFe to the LI2 NisFe. In the 
spin-down subsystem, substantial changes take place in 
the DOS at {E - Ep) > Ep: the part of the d-band ly- 




Fig. 7. Scaled lifetime rkno- x (ekno- — Ef)'^ (upper panel) and 
renormalization factor Zi^na (lower panel) as functions of the 
exciting energy ekncr — Ep for spin-up (light gray points) and 
spin-down (dark gray points) quasiparticles in the Llo NiFe 
ferromagnetic compound. Solid and dashed lines show the cor- 
responding scaled lifetimes averaged over momentum k for a 
given exciting energy. 






N 




Fig. 8. Same as in Fig. [7] but in the LI2 NiaFe ferromagnetic 
compound. 



ing right above Ep has narrowed and the corresponding 
DOS has become higher. This DOS is characterized by the 
comparable contributions of Fe and Ni. At that, under Ep 
spin-down nickel states dominate. As a consequence, the 
momentum-averaged factor Z^ni has changed its behavior 
around Ep (approximately ±1.5 eV), where now there is 
a wide minimum. 
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In Fig. ini we show all our results obtained for the Ni- 
Fe system. As is seen from the figure (see upper panel), 
like in the Co-Fe system the lifetime of spin-up electrons 
in the compounds is longer than that in the bcc Fe and 
the fee Ni. As to the lifetime of spin-down electrons in the 
compounds, it demonstrates a cross between Fe and Ni. 
At that, due to the aforementioned dominant contribution 
of iron states to the DOS in the spin-down subsystem of 
the Llo NiFe, Tknj. in the latter is closer to that in Fe. 
High atomic content of nickel in the LI2 NiaFe results 
in the lifetime Tknj,, whcih tends to be closer to that in 
the fee nickel. Note that in the latter spin-up and spin- 
down holes have practically the same lifetimes (except for 
a small vicinity of Ep). At exciting energies less than —1 
eV, the corresponding lifetimes of holes in the compounds 
are quite similar. 

In order to analyze the spin asymmetry, in Fig.[S](lower 
panel) we also represent the ratio T-f/r^. Again, as in the 
case of the Co-Fe system, the shown data correlate with 
the DOS shape and spin polarization in the vicinity of 
the Fermi level. As is follows from the figure, the largest 
ratio corresponds to the fee nickel, which has the highest 
density of spin-down states at Ep. The ratio in Ni de- 
creases from ~ 13 to ~ 8 away from the Fermi level to 
{E — Ep) « 0.6 eV and becomes smaller than t^/t^ in 
the compounds after (E — Ep) w 1.3 eV. On average, for 
electrons the LI2 NiaFe demonstrates the ratio about 5 
against 4 in the case of the Llg NiFe. The inset of the 
lower panel of Fig. [S] shows the inverse lifetime Tj_/r-|-. The 
presented curves reflect the situation with the lifetimes 
of holes, which can be considered as a gradual transition 
from the bcc iron to the fee nickel. 

Now, we analyze the IMFP of quasiparticles in the fer- 
romagnetic materials considered in this subsection. Fig.lTUl 
shows the obtained results on the IMFP of electrons as a 
function of exciting energy in comparison with the ex- 
perimental data on the attenuation length taken from 
[3], where the spin-dependent transport of electrons in 
NigiFcig films on the GaAs(OOl) surface has been stud- 
ied. As is clearly seen, among the materials considered 
here the compound NisFe with the LI2 structure is char- 
acterized by the largest IMFP of spin-up electrons. As to 
the IMFP of spin-down electrons in this compound, it ap- 
pears to be the smallest at {E — Ep) < 1.3 eV and very 
close to the momentum-averaged A^~| in the bcc Fe and 
the Llo NiFe at {E - Ep) > 2.0 eV. Note that in the 
fee Ni and NisFe with LI2 structure the IMFP of spin- 
down electrons is almost independent of exciting energy. 
Such a behavior of the IMFP as a function of energy is 
in agreement with the experimental data. Moreover, the 
calculated values are quite close to the experimental ones. 
However, our G'^W'^ calculations overestimate the atten- 
uation length of spin-up electrons. 

As well as in the case of the Co-Fe system, there is 
a big spin- filtering effect that is not so strong as, e.g., in 
the DO3 C03FC but nevertheless disagrees with the experi- 
mental observation. As before, we approximate the inverse 
value of the attenuation length of spin-up electrons in the 
LI2 NigFc as 1/A|"^ 4- E^/v^ but with F^ ^ 25 meV. The 



-^ NiFe (Llg) 
-^ NijFe (LI,) 

■i Ni(fcc) 

-i-Fe(bcc) 




-1 1 
E-E^(eV) 

Fig. 9. Scaled lifetime averaged over momentum k for a given 
exciting energy (upper panel) and ratio r-|-/rj_ between the aver- 
aged hfetimes of quasiparticles with different spin orientations 
(lower panel) as functions of the exciting energy for the fer- 
romagnetic compounds NiFe and NiaFe with the Llo and LI2 
structure, respectively, as well as for the bcc Fe and the fee Ni. 
Inset in the lower panel shows the inverse ratio Tj_/T-f for holes. 



^NiFe(Ly 

^NijFeCLy 

i- Ni(fcc) 

i--Fe(bcc) 

NigjFe^g (experiment) 
l/x""'"'' + 25meV 




2.0 2.5 3.0 

£-£^(eV) 

Fig. 10. Dependence of the momentum-averaged IMFP A|J,^^ 
on the exciting energy. The experimental data on the attenu- 
ation length for spin-up and spin-down electrons in NisiFeig 
are taken from [3] 



obtained results are presented in Fig. [TU] We thus have 
achieved an agreement with the experimental attenuation 
length decreasing monotonically with increasing exciting 
energy. This means that some contribution that depends 
weakly on energy (similar to the electron-phonon broad- 
ening) should be taken into account. 



4 Conclusions 

In conclusion, we have presented the spin-polarized G^VF^ 
calculations of the quasiparticle lifetime and the quasipar- 
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tide mean free path caused by inelastic electron-electron 
scattering in ferromagnetic pure metals and compounds 
of the Co-Fe and Ni-Fe systems. Among the compounds, 
we have considered CoFe and NiFe with the B2 and LIq 
structure, respectively. To be closer to the Co- and Ni-rich 
ferromagnetic alloys used in practice, we have also studied 
the DO3 CoaFe and the LI2 NiaFe compounds. Such a set 
of ferromagnetic materials has allowed us to demonstrate 
effects that structure and stoichiometry changes have on 
the quasiparticle properties under study. 

We have found the significant difference between the 
lifetime of electrons with spin up and spin down, which 
is caused by the band-structure characteristic feature of 
the considered ferromagnetic metals and compounds. In 
the DOS, this feature manifests itself as the difference be- 
tween the densities of spin-up and spin-down states in the 
vicinity of the Fermi energy due to the exchange splitting 
of the energy bands. We have represented the resulting 
spin asymmetry of the lifetime by the ratio tj-Ztj, of the 
lifetimes of spin-up and spin-down electrons. We have re- 
vealed that in the compounds the ratio can be noticeably 
larger than that in their constituents. At that, owing to the 
originally large ratio for electrons in the bcc Co as com- 
pared with the fee Ni, the compounds of the Co-Fe system 
demonstrate t^Jt^ that on average is larger than the ratio 
for electrons in the compounds of the Ni-Fe system. On 
the whole, with respect to the lifetime spin asymmetry the 
Co- and Ni-rich compounds are preferable. 

In order to estimate the contribution of the inelastic 
electron-electron scattering to the spin-filtering effect ex- 
perimentally observed in spin-dependent transport of elec- 
trons in the systems under study, we have analyzed the in- 
elastic mean free path of spin-up and spin-down electrons 
as a function of exciting energy. We have shown that the 
mean free path expressed as the product of the lifetime 
and the quasiparticle velocity inherits the lifetime spin 
asymmetry increased by the noticeable difference between 
spin-up- and spin-down-electron velocities. As it follows 
from our calculations, practically within the whole pos- 
itive exciting-energy range considered in the paper the 
compounds of the Co-Fe system possess the longest inelas- 
tic mean free path of spin-up electrons in comparison with 
the compounds of the Ni-Fe system. At that, the inelastic 
mean free paths of spin-down electrons are quite similar 
for both systems. However, this finding disagrees with the 
available experimental observations. The reason underly- 
ing such a discrepancy can relate to decaying mechanisms 
additional to the inelastic electron-electron scattering. Ac- 
tually, both in the Co-Fe system and in the Ni-Fe one, we 
have unambiguously shown that in order to reach a sat- 
isfactory agreement with the available experimental data 
on the attenuation length a contribution to the quasiparti- 
cle decay rate, which depends weakly on exciting energy, 
should be taken into account in addition to the contri- 
bution coming from the inelastic electron-electron scat- 
tering. We believe that the electron-phonon broadening 
contributes significantly to the attenuation length. In or- 
der to clear up this point, spin-polarized calculations of 



the electron-phonon broadening in the considered ferro- 
magnetic compounds should be performed. 
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